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Abstract 

We have mesured the carrier recombination dynamics in InGaN/GaN multiple quantum wells 
over an unprecedented range in intensity and time by means of time-resolved photoluminescence 
spectroscopy. We find that, at times shorter than 30 ns, they follow an exponential form, and 
a power law at times longer than 1 /is. To explain these biphasic dynamics, we propose a simple 
three-level model where a charge-separated state interplays with the radiative state through charge 
transfer following a tunneling mechanism. We show how the distribution of distances in charge- 
separated states controls the dynamics at long time. Our results imply that charge recombination 
happens on nearly-isolated clusters of localization centers. 

PACS numbers: 78.47.-p, 78.47.Cd, 78.55. Cr, 78.67.De, 
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I. INTRODUCTION 



In^,Gai_^.N/GaN heterostructures are essential for the production of optoelectronic de- 
vices, such as blue/ultraviolet lasers and light-emitting diodes, with external quantum effi- 
ciencies approaching 75%." The physical mechanisms that result in such high efficiencies 
are still poorly understood, as high threading dislocation densities due to epitaxial growth 
on mismatched sapphire substrates are belived to induce numerous non-radiative carrier re- 
combination centers in these heterostructures.-"- It is generally assumed that, because of the 
large electron and hole effective masses in InGaN, nanoscopic potential fluctuations induce a 
strong localization of the carrier wavefunction.-'^ Evidence for carrier localization in In-rich 
clusters^ii^ and antilocalization near threading dislocations^'ii'i^ has been reported. Such 
a localization has a profound impact on carrier dynamics, which can be unraveled through 
time- resolved photoluminescence (PL) measurements. The PL decay from InGaN/GaN 
heterostructures is highly non-exponential^>^>^"— but a detailed understanding of these 
particular dynamics is still lacking. 

Early studies concluded that the recombination dynamics could be well described by a 
stretched exponential.-»^>^"— Stretched exponentials are often observed in highly disordered 
materials,— where they are assigned to the summation of monoexponential decays originat- 
ing from individual localization centers with different exciton lifetimes.— The quantum- 
conflned Stark effet has also been invoked to explain the carrier dynamics. It is induced 
in InGaN/GaN heterostructures by the large electric fields that result from the mechanical 
constraints inherently present in epitaxially grown nitrides, combined with high piezoelectric 
coefficients. It leads to non-exponential recombination dynamics, redshift of the emission 
energy with time, and excitation density dependence on the emission energy.— The quantum- 
confined Stark effet is mostly observed in wide InGaN quantum wells (QW), since in QWs 
whose thickness is less than the exciton Bohr radius the e-h overlap integral is not modified 
by the presence of electric fields.— The localization of electrons and holes on separate sites 
has also been reported.— i2ii2I The amphoteric character of the localization centers has led 
to proposals of a bidimensional donor-acceptor like (2D-DAP) recombination model where 
electrons and holes are localized on the opposite sides of the QWs. The 2D-DAP model 
was found to reproduce the PL decay dynamics of InGaN/GaN QWs over three decades 
in intensity.— Recently, the non-exponential PL dynamics of ensembles of localized centers 
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have been linked to the PL intermittency, also called blinking, of individual centers whose 
emission is observed under continuous excitation conditions.— Blinking is observed in a 
variety of systems such as colloidal and self-assembled quantum dots (QDs)^i2^""— and Si 
nanocristals,^^ and is attributed to the formation of metastable, charge-separated (CS) dark 
states^iSa 

In this paper, we present ensemble measurements of the carrier recombination dynam- 
ics in InGaN/GaN QWs over a time window that allows us to follow the PL decay over 
six decades in intensity. We observe a previously unreported slowing down of the decay 
dynamics at times longer than 200 ns that cannot be explained by stretched exponential, 
quantum-confined Stark effet or 2D-DAP models. We propose a three-level recombina- 
tion model involving a dark CS state to explain our results. We conclude that, in thin 
InGaN/GaN QWs, radiative recombination happens on nearly-isolated islands of agglomer- 
ated localization centers. 

II. EXPERIMENTAL DETAILS 

The InGaN/GaN samples were grown on c-plane sapphire substrates by metalorganic 
chemical vapor deposition. Five-period Inj.Gai_j.N/GaN QWs were grown on 2 fim n-GaN 
with a nominal In composition x = 0.2. The thickness of the GaN barriers was 14.0 and 
that of the InGaN wells between 2 and 4 nm. On top of the multiple QW structure, 14 nm 
of p-AlGaN and 100 nm of p-GaN were deposited. 

Time-resolved PL measurements were performed using a Ti:sapphire laser system with 
a repetition rate of 1 kHz and a pulse duration of 40 fs. The second harmonic (390 nm) of 
the laser beam, generated with a type-I /3-BBO crystal, was focussed on the samples with 
an excitation density of about 3 x lO^Wcm"^, corresponding to 6 x 10^^ photons cm^^ per 
pulse. This excitation density is well within the linear regime of our samples. The samples 
were cooled to 9 K in a cold-finger liquid-helium continuous-flow cryostat. As the excitation 
photon energy (3.18 eV) is lower than the low-temperature GaN band gap (3.5 eV), the 
carriers were photogenerated directly in the QWs. The PL signal was analyzed with a 0.3-m 
grating spectrometer and detected with a gated, intensified charge-coupled device. The time 
evolution of the emission spectrum was acquired with a temporal resolution of 5 ns. 
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III. RESULTS AND DISCUSSION 



The PL spectra of our samples showed multimodal emission spanning the range 2.2-3.1 
eV. All samples exhibited non-exponential decay dynamics with similar temporal dynamics. 
The time-resolved PL spectra of the sample with 2 nm QWs are presented in Fig. [H Three 
bands can be distinguished. Band A, at 2.81 eV, has a lifetime shorter than 4 ns. We 
attribute it to QW free-exciton recombination. Band B peaks at 2.74 eV. After a rapid 
energy shift for t < 4 ns, its shape and energy position remain constant with time. At times 
longer than 1 fis, band C can be observed at 2.65 eV. The time-integrated PL spectrum is 
close to that of the time-resolved emission at 10 ns. As band B dominates the emission, we 
focus on its dynamics in what follows. 

Figure[2] shows the decay of the emission from band B integrated from 2.73 to 2.76 eV. 
For times shorter than 30 ns, the intensity decreases rapidly and its dynamics are compatible 
with a stretched exponential. At longer times however, the dynamics slow down and tend to 
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Figure 1. Time-resolved PL spectra of the sample with 2 nm QWs. TIPL is the time-integrated 
emission. Eexc = 3.18eV. 
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Figure 2. PL decay of band B. a) semilogarithmic plot; b) double logarithmic plot. The solid lines 
are the best fits using the stretched exponential (SE), 2D-DAP, and CS models. 



a power law t~" with a = 1.3 for t > 200 ns, as can be seen on Fig. 2(b) The long-lived tail 
on the decay curve, which accounts for about 20% of the time integrated intensity, imposes 
stringent constraints on any dynamical model. 

The absence of spectral shift with time, the invariance of the decay curve with excitation 
density, together with the small thickness of the InGaN QWs, rules out the presence of a 
quantum-confined Stark effect in our samples. Further, the significant proportion of photons 
emited in the long-lived tail implies that the power law it is not caused by a rare occurence 
of localization centers or impurities. In addition, the shape preservation of bands with 
time and the energy independence of the decay dynamics, indicates the presence of a single 
recombination mechanism. 

Within the stretched exponential model, the PL intensity is assumed to follow 



/(t) = /oexp|-(rt)^ 



(1) 



As can be seen in Fig. 2(a), Eq. [T] reproduces our data with the parameters listed in Table 



[J only for t < 30 ns, that is for less than two decades in intensity. Thus, in our samples, the 
stretched exponential model is at best a phenomenological description at short times. 



Table I. Parameters obtained by fitting the stretched exponential (SE), 2D-DAP ans CS models 
to the experimental data. The numbers in parenthesis are the uncertainties on the values of the 
parameters. 
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Values obtained with v = 2 x IQ^^ s ^. 



The 2D-DAP model implies a decay dynamics that follow^ 

poo 

I{t) =2Txn / iy(r)e-^(^)Vdr 
Jo 

X exp <|27m ^~(e-^('-)*-l)rdr|, (2) 

where W{r) = Pe"^^/"^ represents the recombination rate for an electron-hole pair separated 
by a distance r, a is the Bohr radius, and n is the majority carrier saturation density. Eq. 
|2] implies that the temporal shape of the PL decay remains the same in all samples where 
7 = na^ is constant. Such a time-scale invariance was indeed reported by Morel et alr^ 

Figure [2] shows the best fit of our data with Eq. |5J A good agreement is achieved for times 
up to about 200 ns with the optimized parameters listed in Table [H The value obtained for 
7 in our samples is close to that obtained by Morel et al— This indicates that our sample 
has a similar short-time decay dynamics. However, the 2D-DAP model fails to reproduce 
the decay slowing down observed at t > 200 ns. 

Sher et al. introduced a three-level CS model to describe both PL blinking and non- 
exponential recombination dynamics in colloidal quantum dots.^^ It is based on the interplay 
between a bright state (1) that corresponds to an electron and a hole on the same localiza- 
tion center, and a dark state (2) where the carriers are spatially separated, as schematized 
in Fig. 131 Luminescence comes from transitions 1 — > with a rate constant Prad, while tran- 
sitions 1 — > 2 occur with a rate constant Pgsc- Within the CS model, the charge separation 
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Figure 3. CS recombination model. State 1 corresponds to an electron and a hole on the same 
localization center. State 2 corresponds to carriers localized on different centers. Radiative transi- 
tions from state 1 to the ground state occur with rate constant F^ad and charge separation from 
state 1 to state 2 occurs with a rate constant Fesc- 

induces a local modification of the potential landscape that affects the back-transfer 2 — )■ 1 . 
If one assumes that the back-transfer occurs by tunnelling, the distribution of time spent in 
the dark state is 

POO 

R(t) = / f{r)k{r)e-^^'^'dr, (3) 
Jo 

where r is the distance between the carriers, k{r) = ue'"^^ is the tunnelling rate, and /(r) 
is the spatial distribution of charge separation.^* Within the CS model, the PL intensity is 
given by 

m it), (4a) 

where xi is the population of state 1 and 

da;i(t) 



dt 



— i^rad + Tesc) ^lit) 

+ [ TescXiit')Rit~t')dt'. 
Jo 



(4b) 



In the Laplace domain, 



Xl(s) = (s + Trad + Lescfl " -R(s)] 



(5a) 



with 




(5b) 



If an exponential spatial distribution is assumed, /(r) = ee and 




(6) 



where x = rjr and fi = e/rj. R{s) can be evaluated through numerical integration of Eq. [6] 



and I{t) by an inverse Laplace transform of Eq. [5a] using the Gaver-Stehfest algorithm.— 
The CS model predicts biphasic dynamics. At short times, the PL decays exponentially, 
I{t) oc e~^'~°-'^*, while it follows a power law at long times, I{t) ~ t~^^'^^\ Equations I4bl and 
E] indicate that u and T^sc cannot be determined independently. For a given tunnel attempt 
frequency u, Fg^c controls the strength of the long-time power law component with respect 
to the short-time exponential one. 

As shown in Fig. [2l the CS model reproduces the experimental results throughout the 
entire range of the measurement with a vibronic-like attempt frequency— and the optimized 
parameters listed in Table HI However, excellent agreement can still be obtained over a wide 
range of attempt frequencies; our simulations show that Tesc ~ ^rad and z/ > 10^° s~^. 

Within our model, the distribution of distances /(r) controls the recombination dynamics. 
We have tested that only rapidly decreasing functions, such as an exponential, can reproduce 
our results. This implies that carriers are confined in nearly isolated islands of grouped 
localization centers. Thus, one can hypothesize that many localization centers are in fact 
present in a single quantum dot. 

The CS model is based on microscopic physical processes similar to the 2D-DAP model, as 
both models involve charge separation. Within the CS model however, the charge separated 
state is dark and can be recycled from the bright state. It is this recycling that causes 
the decay slowing down observed at long times. As indicated above, the 2D-DAP model 
explains nicely the PL decay shape invariance reported in similar InGaN/GaN samples. The 
CS model also obeys a time-scale invariance law. This can be seen from Eq. |5al which can 
be approximated by^ 



Eq. [7] breaks down for /i = 1 but reproduces well the decay dynamics for /i < 1. It follows 
that the temporal shape of the PL decay is invariant if ^rad^esc^~^ remains constant. Thus, 
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within the CS model, the radiative rate of the bright state Trad is hnked to the escape rate 
to the dark state Tesc and to the back-transfer attempt frequency u. Further experimental 
and theoretical work is required to ascertain the physical origin of this time-scaling law. 

Our results also allow to estimate the energy barrier between the dark and brigth state, 
Eioc, from the back-transfer tunnel rate rj = (8mEioch~'^Y^'^ . Given the composition of our 
wells, one can estimate the electronic effective mass as m = 0.18 mo. The fit of recom- 
bination dynamics gives fi = 0.30. Assuming intra-dot carrier localization, the character- 
istic distance between localization centers is then 1/e ~ 3-6 nm.—^^ One then estimates 
-25-50 meV. 

The proposed model for recombination dynamics in InGaN also seems to be valid in 
a variety of materials such as colloidal QDs^^ or polymer distributed heterostructures.— 
This indicates that the interaction of a charge-separated state with a radiative state plays 
a crucial role in the determination of long-time carriers recombination dynamics in systems 
with some degree of disorder. 

IV. CONCLUSION 

We have measured the time- resolved photoluminescence of InGaN/GaN multiple quan- 
tum wells over an unprecedented range in intensity. The dynamics show a slowing down at 
long times that cannot be explained by currents models. Rather, it can be well modelled 
by a simple three-level system with a dark charge-separated state and a bright radiative 
state. The interplay between the dark and bright states, throught back-transfer tunnelling, 
dominates the recombination dynamics at long times. Our results indicate that the distri- 
bution of distance between charge-localization centers decreases rapidly with distance. This 
implies that the radiative recombination happens on nearly-isolated islands of agglomerated 
localizations centers. 

ACKNOWLEDGMENTS 

This work was supported by the Natural Science and Engineering Research Council of 
Canada (NSERC) and by the Fond Quebecois de la Recherche sur la Nature et les Technolo- 
gies (FQRNT). CS also acknowledges financial support from the Canada Research Chair 

9 



Program. 



* Present address: Universit Europeenne de Bretagne, Prance, INSA, POTON, UMR 6082, P- 

35708 Rennes. 
^ lrichard .leonelli@ umontreal.cal 

^ S. P. Shuji Nakamura and G. Pasol, The blue laser diode (Springer- Verlag, 2000). 

^ B. Gil, Low- dimensional nitride semiconductors (Oxford Univ. Press, 2002). 

^ Y. Narukawa, M. Sano, M. Ichikawa, S. Minato, T. Sakamoto, T. Yamada, and T. Mukai, Jpn. 

J. Appl. Phys. Part 2 46, L963 (2007). 
^ S. J. Rosner, E. C. Carr, M. J. Ludowise, G. Girolami, and H. I. Erikson, Applied Physics 

Letters 70, 420 (1997). 
5 A. Kaneta, M. Punato, and Y. Kawakami, Phys. Rev. B 78, 125317 (2008). 
^ T. Sugahara, H. Sato, M. Hao, Y. Naoi, S. Kurai, S. Tottori, K. Yamashita, K. Nishino, L. T. 

Romano, and S. Sakai, Jpn. J. Appl. Phys. 37, L398 (1998). 

H. Sch5mig, S. Halm, A. Porchel, G. Bacher, J. Off, and P. Scholz, Phys. Rev. Lett. 92, 106802 
(2004). 

^ D. M. Graham, A. Soltani-Vala, P. Dawson, M. J. Godfrey, T. M. Smeeton, J. S. Barnard, M. J. 

Kappers, C. J. Humphreys, and E. J. Thrush, J. Appl. Phys. 97, 103508 (2005). 
^ I. L. Krestnikov, N. N. Ledentsov, A. Hoffmann, D. Bimberg, A. V. Sakharov, W. V. Lundin, 

A. P. Tsatsul'nikov, A. S. Usikov, Z. I. Alferov, Y. G. Musikhin, and D. Gerthsen, Phys. Rev. 

B 66, 155310 (2002). 

^ T. Bartel, M. Dworzak, M. Strassburg, A. Hoffmann, A. Strittmatter, and D. Bimberg, Applied 

Physics Letters 85, 1946 (2004). 
^ A. Hangleiter, C. Netzel, D. Puhrmann, P. Hitzel, L. Hoffmann, H. Bremers, U. Rossow, G. Ade, 

and P. Hinze, Philosophical Magazine 87, 2041 (2007). 
^ Y. Narukawa, Y. Kawakami, M. Punato, S. Pujita, S. Pujita, and S. Nakamura, Appl. Phys. 

Lett. 70, 981 (1997). 

^ Y. Sun, Y.-H. Cho, E.-K. Suh, H. J. Lee, R. J. Choi, and Y. B. Hahn, Applied Physics Letters 
84, 49 (2004). 

10 



M. Pophristic, F. H. Long, C. Tran, I. T. Ferguson, and R. F. Karlicek, Jr., Journal of Applied 
Physics 86, 1114 (1999). 

S. Chichibu, T. Onuma, T. Sota, S. P. DenBaars, S. Nakamura, T. Kitamura, Y. Ishida, and 

H. Okumura, Journal of Applied Physics 93, 2051 (2003). 

T. Onuma, A. Chakraborty, B. A. Haskell, S. Keller, S. P. DenBaars, J. S. Speck, S. Nakamura, 
U. K. Mishra, T. Sota, and S. F. Chichibu, Applied Physics Letters 86, 151918 (2005). 
1^ S.-Y. Kwon, H. J. Kim, E. Yoon, Y. Jang, K.-J. Yee, D. Lee, S.-H. Park, D.-Y. Park, H. Cheong, 
F. Rol, and L. S. Dang, Journal of Applied Physics 103, 063509 (2008). 

P. Lcfcbvrc, S. Kalliakos, T. Bretagnon, P. Valvin, T. Taliercio, B. Gil, N. Grandjean, and 
J. Massies, Phys. Rev. B 69, 035307 (2004). 

H. Scher, M. Shlesinger, and J. Bender, Phys. Today 26, 24 (1991). 

20 P. Bigenwald, P. Lefebvre, T. bretagon, and B. Gil, phys. stat. sol. (a) 216, 371 (1999). 

21 A. Morel, P. Lefebvre, S. Kalliakos, T. Taliercio, T. Bretagnon, and B. Gil, Phys. Rev. B 68, 

045331 (2003). 

22 S. Kalliakos, X. B. Zhang, T. Taliercio, P. Lefebvre, B. Gil, N. Grandjean, B. Damilano, and 
J. Massies, Applied Physics Letters 80, 428 (2002). 

23 P. H. Sher, J. M. Smith, P. A. Dalgarno, R. J. Warburton, X. Chen, P. J. Dobson, S. M. Daniels, 
N. L. Pickett, and P. O'Brien, Appl. Phys. Lett. 92, 101111 (2008). 

24 K. Dunn, J. Derr, T. Johnston, M. Chaker, and F. Rosei, Phys. Rev. B 80, 035330 (2009). 

25 M. Nirmal, B. O. Dabbousi, M. G. Bawendi, J. J. Macklin, J. K. Trautman, T. D. Harris, and 
L. E. Brus, Nature 383, 802 (1996). 

26 M. Kuno, D. P. Fromm, A. Gallagher, D. J. Nesbitt, O. I. Micic, and A. J. Nozik, Nano Letters 

I, 557 (2001). 

2^^ D. Bertram, M. C. Hanna, and A. J. Nozik, Applied Physics Letters 74, 2666 (1999). 

28 M. Tachiya and K. Seki, Appl. Phys. Lett. 94, 081104 (2009). 

29 J. Zhao, G. Nair, B. R. Fisher, and M. G. Bawendi, Phys. Rev. Lett. 104, 157403 (2010). 
3° H. Stehfest, Commun. ACM 13, 47 (1970). 

31 A. C. Morteani, A. S. Dhoot, J.-S. Kim, C. Silva, N. C. Greenham, C. Murphy, E. Moons, 
S. Cina, J. H. Burroughes, and R. H. Friend, Advanced Materials 15, 1708 (2003). 

32 P. Rinke, M. Winkelnkemper, A. Qteish, D. Bimberg, J. Neugebauer, and M. SchefHer, Phys. 
Rev. B 77, 075202 (2008). 

11 



S. Chichibu, T. Azuhata, T. Sota, and S. Nakamura, Appl. Phys. Lett. 69, 4188 (1996). 

1. P. Soshnikov, V. V. Lundin, A. S. Usikov, I. P. Kalmykova, N. N. Ledentsov, A. Rosenauer, 

B. Neubauer, and D. Gerthsen, Semiconductors 34, 621 (2000). 

S. Gelinas, C.-N. Brosseau, O. Pare-Labrosse, I. A. Howard, K. R. Kirov, R. H. Friend, 
R. Leonelli, and C. Silva, unpublished. 



12 



